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We investigate the non-Fermi-liquid (NFL) behavior of the impurity Anderson model (IAM) 
with non-Kramers doublet ground state of the f 2 configuration under the tetragonal crystalline 
electric field (CEF). The low energy spectrum is explained by a combination of the NFL and 
the local-Fermi-liquid parts which are independent with each other. The NFL part of the spec- 
trum has the same form to that of two-channcl-Kondo model (TCKM). We have a parameter 
range that the IAM shows the — In T divergence of the magnetic susceptibility together with 
the positive magneto resistance. We point out a possibility that the anomalous properties of 
U a; Thi_ :E Ru2Si2 including the decreasing resistivity with decreasing temperature can be ex- 
plained by the NFL scenario of the TCKM type. We also investigate an effect of the lowering 
of the crystal symmetry. It breaks the NFL behavior at around the temperature, 8/ '10, where 8 
is the orthorhombic CEF splitting. The NFL behavior is still expected above the temperature, 
(5/10. 
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§1. Introduction 

The material, U x Thi_ x Ru2Si2, seems to show the non- Fermi-liquid (NFL) behavior even in the 
dilute uranium limit flii Its magnetic susceptibility and y-coefficient of the specific heat show the 
— InT divergence at the low temperatures. At the same time the electrical resistivity decreases 
with decreasing temperature. Particular interests have been aroused in the NFL behavior driven 
by the single site effect.^ This material has the tetragonal type structure, and the valence of the 
uranium ion is expected to be mainly, U 4+ (5f 2 ). The crystalline electric field (CEF) ground state 
is also expected to be the non-Kramers doublet state from the measurement of the magnetization. 

Many theoretical studies have been done for the NFL behavior of the dilute uranium compounds. 
Cox has pointed out the possibility that the electronic state of ion with f 2 configuration in a specific 
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situation can be mapped on the S= 1/2- two-channel Kondo model (5=1/ 2-TCKM)Jl!) which i 
known to show the NFL behavior. The required situation is that the lowest CEF state 

of the f 2 configuration is the non-Kramers doublet state, and two kinds of conduction electrons, 
both of which have doubly degenerate time reversal pair, can hybridize with the f state. However, 
there are two remained problems in the application of the £=l/2-TCKM to the NFL behavior 
of U x .Thi_ x Ru2Si2: First, the S=l/2-TCKM can not explain the decrease of the resistivity with 
decreasing temperature. Usually opposite temperature dependence is expected when the coefficient 
of the — In T term of the magnetic susceptibility is sizable^© Secondly, the mapping is done 
assuming very restricted situation. It is not clear whether the NFL state is stable or not when the 
realistic situation of the CEF states is taken into account. For the first problem we have modified 
the TCKM, and have proposed the extended two-channel Anderson model (ETCAM) that shows 
similar behaviors of the resistivity and the magnetic susceptibility to those of the experiment 
However, applicability of the ETCAM to the uranium problem is not obvious at present, because 
the intimate one-to-one mapping of electronic state has not been verified. In this paper, we check 
whether the anomalous properties of U x Thi_ a; Ru2Si2 can be explained or not based on a more 
realistic model. 

For the second problem Koga and Shiba have studied the stability of the NFL state by taking 
into account the first excited state of the f 2 configurationJli'0llZl' They applied the numerical 
renormalization group (NRG) method, and has shown that the NFL state is still stable, if the 
CEF splitting is large. In previous paper, we have also studied the stability of the NFL state 
based directly on the impurity Anderson model (IAM),® because preceding studies have been 
done applying the Schrieffer- Wolff transformation by assuming the strong correlation limit. It has 
been shown that the ground state property can change to the NFL state from the usual Kondo 
type singlet state when the intensity of the hybridization is gradually weakened, even though all 
CEF states are taken into account, and the valence is apart from f 2 . But details of the low energy 
properties of the NFL state have not been clarified. 

The first purpose of this paper is to analyze the low energy spectrum of the IAM given by 
the NRG method in detail.© We examine the model that the single electron orbitals split into 
three doublets under the tetragonal CEF. The CEF ground state is the non-Kramers doublet of f 2 
configuration. In §3 we show that the low energy spectrum is explained by a combination of two 
components which are given by the NFL and the local- Fermi-liquid (LFL) fixed point Hamiltonians, 
and which are independent with each other. The NFL part of the low energy spectrum has the same 
form to that of the TCKM. The second purpose is to consider whether the anomalous behavior in 
Ua;Thi_ x Ru2Si2 is explained from the IAM or not. In §4 and 5 we show that the large — InT term 
of the magnetic susceptibility can be expected together with the positive magneto resistance at very 
low temperature, when the hybridization is not so weak. The positive magneto resistance indicates 
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that the temperature dependence of the resistivity decreases with decreasing temperature. The 
third purpose is to study the effect of the releasing process of the residual entropy which is inherent 
of the S=l/2-TCKM. In §6 we examine the NFL behavior in the magnetization by breaking the 
NFL state with lowering the crystal symmetry. 

§2. Model 

We consider the following IAM under tetragonal CEF, 

H = H f + H cf + H c , (2.1) 
n f = J2 £ ( r ) n /r 7 + J Yl ™/r 7 n /rv 

~~ Iq X! [jr7r'7'Jr"7"r'"7'"/r 7 /r'7'/r» 7 »/r'"7'" 
{r 7 } 
35n/" 
4~ _ 

W C / = EE (Vfl y c kVl + h.c.) , (2.3) 

k F'y 

Hc = mi £ fc4r 7 c fcr7> ( 2 - 4 ) 

k T7 

where /r 7 (c^r-y) is the annihilation operator of f-electron with the 7-th component of T-irreducible 
representation (conduction electron with wave number k), and e(T), U and / denote the single 
f-electron energy, Coulomb and exchange interaction constants, respectively. For simplicity we 
assume the large spin-orbit interaction of f-electron, so only the j = 5/2 orbitals are considered. 
The orbitals split into a quartet, and a doublet, |/ 1 r^ 2 \ ±) = ±a/1/6| ± 5/2) =p-\/5/6| =F 3/2) 
under the cubic CEF where m in \m) of the right hand side of the equation is the magnetic 
quantum number of j& The quartet splits into two doublets again, [f 1 ^ 1 ^, ±) = ±^/5/6| ± 5/2) 
±a/1/6| ^3/2) and |/ 1 rg, ±) = | ±1/2) under the tetragonal CEF. The Coulomb and the exchange 
interactions are assumed to be the j-j coupling type. Usually, the multiplet structure of U ion 
is approximated by the L-S coupling scheme. So the present j-j coupling scheme seems to be 
not applicable to U ion at first glance. However, we note the ground multiplets have same total 
angular momentum for both coupling schemes. We expect qualitative features of the low energy 
states within the ground multiplet will be not so changed. It is assumed the f-electron hybridizes 
with the conduction electron which has the same component of the irreducible representation. The 
quantities, V and denote the hybridization matrix and the band energy, respectively. It is 
also assumed the band is extending in energy from — D to D with constant hybridization matrix, 
r = ttV 2 /2D. The energy unit and the origin of the energy are chosen as D = 1 and the Fermi 
level, respectively. 

We rewrite the Hamiltonian into a form to suite the NRG calculation. First the conduction band 
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is discretized by the logarithmic mesh to give good sampling to states near the Fermi energy. 
Next the Hamiltonian is transformed into an expression given by the shell orbits, 



H = lim H L , (2.5) 

L— >oo 

n L = n f + J2 (^Vf^sory + h.c.) + H° L , (2.6) 

L-l 

«"=EE^ (4+ir 7 ^r 7 + h.c.) , (2.7) 

£=o r 7 

where sgr-y is the annihilation operator of l-th shell state with the 7-th component of the T- 
irreducible representation. The hopping energy, tg, between shell states is given by tg = D(l + 
A _1 )A _ ^/ 2 ^/2, where A(> 1) is the discretization parameter, and £g tend to 1 when i increases. 
The quantity A/y is the correction factor of order 1 for the discretization. We first diagonalize the 
TLf term, and then by adding the shell state successively from I = 0, we diagonalize the series of 
cluster Hamiltonian, {TCl}, recursively. At each step we retain about the 500 lower energy states to 
the next step. This number is not so large, but the obtained eigen states in the low energy region 
seem to have enough accuracy for qualitative discussions. When the only 300 states are retained 
at each step, the low energy eigen states are not changed essentially. 

The valence of U ion in U x Thi_ x Ru2Si2 is not so clear, so we consider two cases: one is the case 
that the valence fluctuation between and Sf 1 configurations becomes dominant, and the other 
is that between bf 2 and 5f 3 configurations. We call the former a f 2 - f 1 dominant fluctuation case, 
and the latter a f 2 - f 3 dominant fluctuation case hereafter. 

§3. Analysis of the Low Energy Fixed Point 

3.1 Flow chart of energy levels of f - / dominant fluctuation case 

We first consider the f 2 - f 1 dominant fluctuation case. We choose the parameters; e(T^) = —0.9, 
e(r 6 ) = -0.75, e(rj 2) ) = -0.5, U = 0.6 and 1 = 8. The lowest CEF state becomes the non- 
Kramers doublet of the f 2 configuration which has mainly the character of J = 4, Ts-irreducible 
representation of D4-group, and the energy of the state is —1.446. The state is given as |/ 2 rs, ±) 
= a\ ± 1/2,t3/2) + b\± 3/2,^5/2) + c| ± 5/2, ±1/2), where a = -0.3395, b = -0.1874 and 
c = 0.9218. We have used the notation, \m,m') = fmfm'\®)' wn ere m and m' are the magnetic 
quantum numbers of j. The excited states from the first to the forth are the singlet of the f 2 
configuration with the energies, -1.442, -1.286, -1.208 and -1.167, respectively. These states have 
mainly the character of J = 4, and are the T4, Ti, T2 and r3-irreducible representations of D4- 
group, respectively. The other configurations, f° and from f 3 to f 6 have higher energy. 

In our previous work we have shown that three types of the ground state; the doublet, the CEF- 
singlet-like and the f°-singlet-like ground states appear successively when the hybridization strength 
is increased by fixing the energy of the CEF statesS The low energy spectrum of the last two 
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follows the LFL theory, however, the first one can not be explained by the LFL theory. In this work 
we concentrate on the weak hybridization case and analyze the spectrum of the NFL state. In Fig. 1 
the lower eigen energies of each cluster Hamiltonian, TCl, for the odd NRG step, L, are shown. The 
energies are renormalized by at each step. The renormalized energy levels of the states which 
has the same charge, Q, and the magnetic quantum number, M, change smoothly, and we call 
the figure the flow chart of the renormalized energy levels (FCEL) hereafter. The quantity, M, is 
defined by using modulo 4 as defined in the caption of Fig. 1 because of the tetragonal symmetry. 
The renormalized energy levels tend to fixed values as L increases beyond a step L ~ 23. The 
right hand side of the figure gives the energy spectrum at the low energy fixed point, and the 
hopping energy, ti-i, at L ~ 23 indicates the measure of the energy scale for the cross over that 
the system goes into the low energy fixed point. Each state is indicated by i (the sequential number 
of states from low to high energy), Q and M. The ground state {i = 1,2) is doublet, and the first 
excited states are two doublet states with one particle excitation (i = 3, 5) and one hole excitation 
(i = 4, 6). The second excited states (i = 7 ~ 10) are the one hole excitation, and they have larger 
energy than twice energy of the first excited states. Therefore these low energy states can not be 
explained by the LFL theory. 

We compare the low energy states with those of 5'=l/2-TCKM which are derived by the conformal 
field theory (CFT)© The low energy states of S'=l/2-TCKM are the NFL states, and the energies 
for primary states are known as conformal towers, 



where Qc,j,jf and 5 P denote the freedoms of the charge, the spin and the flavor-spin, and the 
potential scattering, respectively. As seen from the particle hole symmetry of the first excited 
states obtained by the NRG calculation, 5 P becomes zero. In Table I we compare the low energy 
states from the NRG calculation and those of the CFT. The renormalized energy, E, in the column 
of the NRG is multiplied by a constant factor, r*, so the energy of the first excited states agrees 
with that of the CFT. The correspondence of the charge and the degeneracy between both results 
is perfect. The difference in energy is small for states in lower energy region with r*E ~ 0.5, and 
it increases for ones in larger energy region, r*E ~ 1.0. But the difference decreases when the 
discretization parameter, A in the NRG method, is decreased as shown in Fig. 2. It will disappear 
in the continuum limit, A — > 1. This behavior is the same as that of the S*=l/2-TCKM. It can 
be concluded the NFL part of the low energy spectrum in Fig. 1 has the same form to that of the 
5=l/2-TCKM. 

In the NRG calculation, we have low energy states which are not listed in Table I. They are 
explained as the states accompanied by the extra excitations obeying the LFL theory. In Table II 




(3.1) 



5 



(16, 








-1 


( i 
( 15 


l 




i , 
, 


2M) 
1) 




1.0 


(14, 








3 


( 13 





, 


3) 




- ~ : 


(12, 








1 


(H 





, 


-1) 














(10 


-1 


, 


-4) 


~W 




( 9, 
( 7, 


-2 
-2 






1 

"3 


( 8 
( 6 


-2 
-2 


, 
, 


-1) 
3) 


o; 

J3 - 












( 5 


- 1 


, 


0) 














( 4 


- 1 


, 


0) 














( 3 


-1 


, 


-4) 






( 2, 


- 1 





"2 


( 1 


- 1 


, 


2) 




9.% = 




Q ,2M) 



11 21 
Shell No L(odd) 



0.0, 
0.0, 



19 , 


-1 


, 


-4) 


20 


- 1 


17, 


-1 


0, 


2) 


18 


1 


15 , 


-3 


, 


2) 


16 


1 


13 , 


-3 


, 


-2) 


14 


- 1 


1 1 , 


-1 


0, 


0) 


12 


- 1 


9 , 


-2 


, 






-2 


7, 


-2 


0, 


-I) 


l l 


-2 


5 , 





, 






-2 


3 , 





, 


-11 


I 


-2 


1 , 


-1 


0, 


2)( 2 


- 1 



(24, 0.0, 
(22, 0.0, 

10, 
0, 



?] 
-4) 

2 
-2 
-2 

0) 

!j 

1] 
-2) 



Fig. 1. FCEL of the f 2 - f 1 dominant fluctuation case for the odd renormalization step, L. (In this case the number 
of shell orbits for the conduction electrons is even, L + 1. For example, number of shell orbits is two for L — 1.) 
The energies are renormalized by at each step. The discretization parameter A — 8 is used, and about 500 

states are retained at each step. The parameters of IAM are; e(r^) = —0.9, e(re) = —0.75, e(T^) = —0.5, 
U — 0.6, 1 = 8 and V 2 /2 = 0.035. The occupation number of the f-electron is 1.92. Each state is labeled by the 
index, i (sequential number of states from low to high energy). The charge of each state is denoted by Q, which is 
defined as total electron number minus half of total orbital number, 3(L + 2). The magnetic quantum number, M 
is the index for grouping of total magnetic quantum number, J z . When the quantity, 2J Z changes from —4 + 8n to 
3 + 8n, where n is some integer, the quantity 2M varies from —4 to 3. The renormalized energies of states which 
has the same Q and M are connected as L changed. 
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Fig. 2. A dependence of the energy belongs to the NFL part at the low energy fixed point. The symbols denote 
the renormalized energy, r*E, given by the NRG calculation. The parameters are the same as those in Fig. 1 
except V. The hybridization width is chosen as V 2 /2 = 0.0433, 0.0350 and 0.0310 for the cases of A = 4, 8 and 
16, respectively. The dashed lines shows the energies which are expected by the CFT for the S=l/2-TCKM. The 
factor, r*, is chosen so each energy of the first excited state agrees with 1/8; r* = 0.5787, 0.5924 and 0.619 for 
A = 4, 8 and 16, respectively. 
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Table I. Comparison between the states belong to the NFL part for the odd NRG step of the f 2 - f 1 dominant 
fluctuation case. The states at the low energy fixed point from the NRG calculation and the low energy states for 
the S'=l/2-TCKM expected from the CFT are listed. The parameters for the NRG calculation are the same as 
those in Fig. 1. The states in the column of NRG are the results at L = 31. Index in the column of NRG is the 
same in Fig. 1. The charge Q — Qo, which is defined from in Fig. 1 and Qo = — 1 of the ground state, coincides 
with Qc- The degeneracy of each line in the column of CFT is given by (2j + l)(2j/ + 1), and it coincides with 
the degeneracy given by NRG. The energy, E, is multiplied by a factor r* = 0.5924 which is defined so to agree 
the energy of the first excited state from the NRG and the CFT results. In the column of CFT the ground state 
energy, 3/16, is subtracted. 



NRG CFT 



index 


Q - Qo dc 


gen. 


E 


r*E 


Qc 


3 


jf Ecfi 


'£/ttvf 


1, 2 





2 













1/2 








3, 5 


1 


2 


0.211 




0.125 


1 





1/2 


1/8 


4, 6 


-1 


2 


0.211 




0.125 


-1 





1/2 


1/8 


11, 12, 14, 17, 19, 20 





6 


0.889 ~ 0.891 


0.527 - 


- 0.528 





1/2 


1 


1/2 


13, 15 


-2 


2 


0.890 




0.527 


-2 


1/2 





1/2 


16, 18 


2 


2 


0.890 




0.527 


2 


1/2 





1/2 


33, 34, 39, 40, 43, 44 


1 


6 


1.17 - 1.18 


0.693 - 


< 0.699 


1 


1 


1/2 


5/8 


35-38, 41, 42 


-1 


6 


1.17 ~ 1.18 


0.693 - 


< 0.699 


1 


1 


1/2 


5/8 


55-58 





4 


1.77 




1.05 





3/2 





1 


79, 80 





2 


1.87 




1.11 





1/2 





1 


113, 114, 117-120 


-2 


6 


2.05 - 2.06 


1.21 


~ 1.22 


-2 


1/2 


1 


1 


115, 116 





2 


2.06 




1.22 





1/2 





1 


121-126 





6 


2.06 




1.22 





1/2 


1 


1 


153-158 


2 


6 


2.10 - 2.11 


1.24 ■ 


- 1.25 


2 


1/2 


1 


1 


we show the list of these states. 


The states, i = 7 


~ 10 and i = 


21 ~ 24, 


are 


interpreted 


as the 



single particle excited state of the LFL part. The former (latter) are classified to the one hole 
(electron) excitation of the occupied (unoccupied) fermion orbitals with energy —0.884 (0.988). 
From the analysis of the magnetic quantum number the symmetry of the fermion orbitals of the 
LFL part is determined as Tj. The NFL part for these eigen states is in its ground state, i = 1 and 
2. For the states from % = 25 to 52 in the Table II the NFL part is in the excited states, i = 3 ~ 6. 
Each state can be explained as the product of the NFL state and the LFL state. One can see the 
perfect correspondence of the charge and the magnetic quantum number between the NRG and the 
expected results. The energy in the column of LFL is estimated from the sum of the energies of the 
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NFL state and the LFL part.EJ These energies give good agreement with those of the NRG, and 
thus the NFL and LFL parts seem to be almost independent with each other. The states, % = 53 
and 54, are interpreted as the two hole excited state of LFL part. The energies of the states agree 
with the sum of the energies of two holes as seem from Table II. This is the characteristic features 
of the LFL excitation. We have analyzed about 150 states as the product of the NFL state and 
the LFL one. In these state the particle-hole pair excitation or the two particles excitation of the 
LFL part are includes. It is concluded the low energy spectrum at the fixed point is explained by 
a model that is the combination of the NFL fixed point Hamiltonian and the LFL one which are 
independent with each other. However, we note the two parts have almost the same low energy 
scale because they tend to the low energy fixed point almost simultaneously at about L ~ 23 as 
seen from Fig. 1. We note the spectrum in the NFL part have the particle- hole symmetry, but 
there is the particle-hole asymmetry in the LFL part, so the total spectrum has the asymmetry. 

In Fig. 3 we show the FCEL for the even renormalization step. The parameters are the same 
as those in Fig. 1. The energy flow is very complicated in the intermediate NRG step, but, the 
excitation spectrum at the larger step {L > 26) for the low energy fixed point is simple. The 
spectrum is explained by the same model as that for the odd renormalization step discussed above, 
that is the combination of the NFL fixed point model and the LFL one. The spectrum of the 
NFL part at the low energy fixed point has the same form as that for S'=l/2-TCKM as shown 
in Table III. The energy, E, has the same quantity as that of the odd NRG step. This is the 
characteristic behavior of the NFL state. In Table IV the low energy states including the excitation 
of the LFL part are shown. The states given by the NRG calculation can be explained by the 
combination of the LFL and the NFL states as shown in the column of LFL. The fermion orbitals 
have the Tj symmetry. Only one orbital above the Fermi level is found even though we analyze 
about 150 states. The single hole excitation is expected to have an energy larger than 2.1. 

There are three kinds of orbitals for the conduction electrons, and Tq. Two of them will 

contribute to the NFL state, and remaining one will be left as the LFL part. To find which orbitals 
are related to the NFL state, we consider a fictitious model that the conduction electrons with the 
symmetry are removed from the original model given in eq. (2.1). In Fig. 4 the FCEL of the 
fictitious model is shown. At the low energy fixed point the energy spectrum does not depend on 
whether the NRG step is odd or even. The low energy spectrum is almost the same as that of the 
NFL part of the original model. There are the small particle-hole asymmetry in the spectrum of 
the fictitious model. All the states at the low energy fixed point can be analyzed as the conformal 
towers of the S=l/2-TCKM with the potential scattering term. The conduction electrons with Fq 
and T^ symmetries contribute to the NFL state, and the electrons with to the LFL state. 
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Table II. Comparison between the states belong to the LFL part for odd NRG step of the f 2 - f 1 dominant fluctuation 
case. The states at the low energy fixed point from the NRG calculation and the low energy states which is expected 
by the LFL theory are listed. The parameters for the NRG calculation are the same as those in Fig. 1. The 
quantities, Q, M and energy in the column of NRG is the results for L — 31 in Fig. 1. All the states in NRG can be 
explained by the product of the NFL state and the single particle state which is indicated in the column of LFL. 
The NFL state is indicated by [NFL(i)], where i is the index of the state shown in Fig. 1 and in Table I. The single 
particle state is indicated by F^ and r",^, where the plus (minus) sign means the particle (hole) excitation. The 
two hole excited state is indicated by [r~^] 2 . The excitation energies, E^rij,) and E^T^.), is given as 0.988 and 
0.844, respectively, to agree with those in the column of NRG. The energy in the column of LFL is estimated from 
the sum of the NFL state energy and the single particle excitation energy. 



NRG 

index Q 2M energy 



7, 9 


-2 


-1, 


3 


0.884 


8, 10 


-2 


1, 


3 


0.884 


21, 23 





-1, 


3 


0.988 


22, 24 





1, 


3 


0.988 


25, 29 


-2 


0, 


-2 


1.096 


26, 30 


-3 


o, 


2 


1.096 


27, 31 


-1 


2, 





1.096 


28, 32 


-1 


-2, 





1.096 


45, 49 


1 


0, 


2 


1.200 


46, 50 


1 


o, 


2 


1.200 


47, 51 


-1 


-2, 





1.201 


48, 52 


-1 


2, 





1.201 


53 


-3 




2 


1.768 


54 


-3 




-2 


1.768 



LFL 

state energy 



[NFL(l)]®r- (1) 


0.884 


[NFL(2)]®r- (1) 


0.884 


[NFL(l)]®r+ (1) 


0.988 


[NFL(2)]®r+ (1) 


0.988 


[NFL(4)]®r- (1) 


1.095 


[NFL(6)]®r ni) 


1.095 


[NFL(3)]®r- (1) 


1.095 


[NFL(5)]®r ni) 


1.095 


[NFL(3)]®r+ (1) 


1.199 


[NFL(5)]®r+ (1) 


1.199 


[NFL(4)]®r+ (1) 


1.199 


[NFL(6)]®r+ (1) 


1.199 


[NFL(i)Mr- {1) ] 2 


1.768 


[NFL(2)Mr- {1) ] 2 


1.768 



When we consider a fictitious model that the conduction electrons with the Tq or ly symmetry 
are removed, the low energy states are explained by a Ising type model: One localized Ising spin 
couples with the conduction electrons through the exchange interaction. 



3.2 Flow chart of energy levels of f - f dominant fluctuation case 
For the f 2 - f 3 dominant fluctuation case we choose the parameters; s(T^) = —1.2, s(Tq) = —1.05, 

(2) 

e(T)- ) = —0.8, U = 0.3 and 1 = 8. The f-levels is chosen to be deeper and the Coulomb interaction 
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( i , Q , 2M) 



33,-2.0, 0) (34,-4.0,-2 
31,-4.0, 2)(32,-2.0,-2 
29, 0.0,-81(30,-8.0, 



( i, Q , 2M) 
(IB, -2.0, 2) 




s (27,-2.0, 2) (28,-2.0,-4 
(25,8.0,-4) (26, 0.0, 2 



(17,-2.0,-2) (16,-2.0,-4)\ 
(15, 1 .0,-3) (14, 1.0, 3) \ 

(13,-2.0, 0)s \ 
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Fig. 3. FCEL of the f 2 - f 1 dominant fluctuation case for the even renormalization step L. The parameters for the 
NRG calculation are the same as those in Fig. 1. See the caption of Fig. 1. 

to be smaller than those for the f 2 - f 1 dominant fluctuation case with fixing the CEF splitting of 
f-level and the exchange interaction. The main states of the f 2 configuration are almost the same as 
those of the previous case. The lowest CEF state is the non-Kramers doublet which has the energy, 
-2.346, and it has mainly the character of J = 4, r^-irreducible representation of D4-group. The 
first to the forth excited states are the singlet which has mainly the character of T4, T\, T2 and 
Ts-irreducible representation of D4-group, respectively. The energy of each state is -2.342, -2.186, 
-2.108 and -2.067, respectively. The main states of the f 3 configuration have mainly the character of 
J = 9/2. They are the Tq doublet, the doublet and the another doublet states with energies 
-1.958, -1.927 and -1.907, respectively. When the hybridization matrix T = ttV 2 /2 is changed from 
0.030-7T to 0.0387T for A = 8, the occupation number of the f-electron varied from 2.233 to 2.255. 

In Fig. 5 (a) the FCEL for the odd NRG step is shown. Each energy spectrum at the low energy 
fixed point is explained by the combination of the NFL fixed point Hamiltonian and the LFL one 
which are independent with each other also in this case. In Table V we show the comparison 
between the low energy states in the NFL part of the NRG result and those of the conformal 
tower of S'=l/2-TCKM. The first excited state is the one-particle excitation doublet (i = 3,4), and 
the second excited state is the one-hole excitation doublet (i = 5,6). There exists the particle- 
hole asymmetry, so the multiplying factor, r* = 0.592, and the potential scattering, 5 P = —0.235, 
are determined from the first and the second excited states. The states of the NRG result show 
good correspondence to the result of CFT. We note the difference of the energies has comparable 
magnitude as that found in Table I. It decreases when we use A = 4, and it will be removed in 
the continuum limit, A — > 1. The states accompanied by the LFL excitation are not listed in the 
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Table III. Comparison between the states belong to the NFL part for even NRG step of the f 2 - f 1 dominant 
fluctuation case. The states at the low energy fixed point from the NRG calculation and the low energy states for 
the S'=l/2-TCKM expected from the CFT are listed. The parameters for the NRG calculation are the same as 
those in Fig. 1. The states in the column of NRG are the results at L = 30. The multiplying factor is r* = 0.5924. 
In the column of CFT the ground state energy, 3/16, is subtracted. See the caption of Table I. 
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table. The LFL states are explained as the excitations of the orbital with the 17V symmetry. This 
situation is common to the f 2 - f 1 dominant fluctuation case. 

In Fig. 5 (b) we show the FCEL for the even NRG step. The NFL part of the low energy states 
are listed in Table VI. The quantities, r* and 5 P , don't depend on the even or oddness of the NRG 
step. The LFL part is not listed in the table. The first excited states which are indicated by the 
index, i = 3 ~ 6, in the right hand side of the figure are the one hole excitations of the fermion 
orbital with Tj symmetry from the ground state doublet in the NFL part. The second excited 
states are the two holes excitations from the ground state. As seen from the figure there are the 
one hole excitations and the two holes excitations of the LFL part from each NFL state. 

When we consider the fictitious model that the conduction electrons with symmetry are 
removed, the spectrum at the low energy shows a similar NFL behavior of the original model. In 

(2) 

other cases; the fictitious model removing the conduction electrons with the Tq or T\ symmetry, 

(2) 

we have the Ising type fixed point. Therefore the Tq and T\ components contribute to the NFL 
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Table IV. Comparison between the states belong to the LFL part for even NRG step of the f 2 - f 1 dominant 
fluctuation case. The states at the low energy fixed point from the NRG calculation and the low energy states 
which is expected by the LFL theory are listed. The parameters for the NRG calculation are the same as those in 
Fig. 1. The quantities in the column of NRG are the results for L = 30 in Fig. 3. The NFL state is indicated by 
[NFL(i)], where i is the index of the state shown in Fig. 3 and in Table III. See the caption of Table II. 
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-1 
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part and the to the LFL part. This fact is the same as that of the f 2 - f 1 dominant fluctuation 
case. 

We note when the original model with the large hybridization strength is considered, the low 
energy states are explained by the LFL theory The lowest energy state at each NRG step changes 
from the crystalline field doublet (L = —1) to the singlet state (L > 15) in the FCEL of the odd 
step. The two electrons are bounded at the intermediate step, L = 15, so the singlet state has the 
character similar to that of the f 4 singlet state. This contrasts with the f°-singlet-like ground state 
in the f 2 - f 1 dominant fluctuation case with the large hybridization strength. 

From the analysis of this section we can deduce the following conclusions: The low energy fixed 
point properties of the present model can be given by the combination of the two independent parts, 
one is the NFL component and the another is the LFL component. The two parts have almost the 
same low energy scale, and in addition the separation of these occurs at almost the same energy 
region with the low energy scale. The conduction electrons with Tq and T y 7 symmetry contribute 
to the NFL part and the F^p symmetry to the LFL part. The above facts are seen commonly both 
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Fig. 4. FCEL of the fictitious model for the odd (a) and even (b) renormalization steps L. The parameters for the 
NRG calculation are the same as those in Fig. 1. See the caption of Fig. 1. 



for the f 2 - f 1 and the f 2 - f 3 fluctuation cases. 
§4. Magnetization 

We calculate the temperature dependence of the magnetization of f-electrons, M = (mj) = 
(J2m=-j m frnfm}, by adding the Zeeman term, TL' = —rrifH z , to eq. (2.1), where m is the magnetic 
quantum number of j. The thermal average at the temperature, Tl, is calculated by the eigen states 
of the Lth cluster as, 

M{TL) ~ Trexp { -(H L + HVT L } ' ^ 
T i = i(l + A- 1 )A( i - 1 )/ 2 /A (4.2) 

where (3 ~ 2 is a parameter for suiting the temperature to the eigen energies of each NRG step. 

In Fig. 6 (a) we show the magnetizations which are normalized by the applied magnetic fields. 
These quantities correspond to the magnetic susceptibility when the field is weak enough. The pa- 
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Fig. 5. FCEL of the f 2 - f 3 dominant fluctuation case for the odd (a) and the even (b) renormalization steps L. The 
discretization parameter A = 8 is used, and about 500 states are retained at each step. The parameters of IAM 
are; e(r^. 1) ) = -1.2, e(r 6 ) = -1.05, e(r^ 2) ) = -0.8, U = 0.3, I = 8 and V 2 /2 = 0.036. The occupation number of 
the f-electron becomes 2.21. See the caption of Fig. 1. 



rameters of the IAM are chosen so that the valence fluctuation between the f 2 and f 3 configurations 
becomes dominant. We use A = 8 and about 400 states are kept at each NRG step. Even when 
we keep about 500 states, the results of the magnetization are not changed essentially. The nor- 
malized magnetizations have almost same value, when the temperature is higher than the Zeeman 
energy. For very weak magnetic field there is a temperature region that the magnetization follows 
the — InT dependence. When the temperature decreases, the magnetization saturates. However, 
— In T temperature region extends with decreasing the applied magnetic field. This behavior can 
be seen clearly from Fig. 6 (b) which shows the differential of the magnetizations by the logarithm 
of the temperature. The magnetic susceptibility follows — In T dependence at low temperatures. 
This temperature dependence is the characteristic behavior of the NFL state for S'=l/2-TCKM 
type. The inverse of the coefficient of — InT term is estimated to be about 10~ 7 . This means that 
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Table V. Comparison between the states belong to the NFL part for the odd NRG step of the f 2 - f 3 dominant 
fluctuation case. The states at the low energy fixed point from the NRG calculation and the low energy states for 
the S'=l/2-TCKM expected from the CFT are listed. The parameters for the NRG calculation are the same as 
those in Fig. 5. The states in the column of NRG are the results at L = 31. The multiplying factor, r* = 0.5920, 
and the potential scattering, S p — —0.2350, are defined to agree the energies of the first and the second excited 
states from the NRG and the CFT results. In the column of CFT the ground state energy, (tt5 p /2) 2 + 3/16, is 
subtracted. See the caption of Table I. 









NRG 










CFT 




index 


Q-Qo 


degen. 




E 


r*E 


Qc 


3 


3! 




1, 2 





2 













1/2 








3,4 


1 


2 




0.148 


0.0876 


1 





1/2 


0.0876 


5, 6 


-1 


2 




0.274 


0.162 


-1 





1/2 


0.162 


7, 8 


2 


2 




0.749 


0.443 


2 


1/2 





0.425 


9-14 





6 


0. 


889 ~ 0.892 


0.524 ~ 0.528 





1/2 


1 


1/2 


23, 24 


-2 


2 




1.03 


0.610 


-2 


1/2 





0.575 


29-34 


1 


6 




1.09 ~ 1.10 


0.645 - 0.648 


1 


1 


1/2 


0.588 


47-52 


-1 


6 




1.26 ~ 1.27 


0.748 ~ 0.749 


-1 


1 


1/2 


0.662 



the energy scale which characterizes the strength of the — In T term of the susceptibility is about 
ICT 7 , and is very small. Therefore the present model gives sizable — In T term. 

The saturation of the magnetization relates to a break of the NFL state under the magnetic field 
as follows: In Fig. 7 the FCEL for the odd NRG step under the magnetic field, H z = 1.0 x 1CP 9 , is 
shown. The energy flow in the region below L = 15 is almost identical to that in Fig. 5(a) which 
is the FCEL without the magnetic field. The spectrum for L = 15 ~ 21 splits gradually due to the 
magnetic field. However, the behavior is very similar to that of the NFL spectrum at the low energy 
fixed point in Fig. 5 (a). The hopping energy at L = 21, tL-i ~ A~ 10 ~ 10 -9 , has the same order 
to that of the Zeeman energy. We note that the magnetization shows the saturation deviated from 
the — InT behavior at about T ~ 10~ 10 for the case of H z = 10 -9 . The flow in the region above 
L = 21 goes to the LFL spectrum at the low energy fixed point. The first and second excited states, 
i=2 and 3, where i is the index in the figure, are the one electron excitations, respectively. The 
third (i=4) and fourth (i=5) excited states are the one hole excitations. The excited state, i=6, is 
the two electrons excitation which is explained by the combination of the excitations for the states, 
i=2 and 3. The excited states, i = 7 ~ 9, are explained by the electron-hole pair excitations, and 
other high energy states are also explained by the LFL fixed point model. The NFL state breaks 
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Table VI. Comparison between the states belong to the NFL part for the even NRG step of the f 2 - f 3 dominant 
fluctuation case. The states at the low energy fixed point from the NRG calculation and the low energy states 
for the S'=l/2-TCKM expected from the CFT are listed. The parameters for the NRG calculation are the same 
as those in Fig. 5. The states in the column of NRG are the results at L = 30. The multiplying factor and the 
potential scattering are: r* = 0.5799 and S p = —0.2393, respectively. In the column of CFT the ground state 
energy, (tv8 p /2) 2 + 3/16, is subtracted. See the caption of Table I. 
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due to the magnetic field, and the low energy states can be explained by the LFL theory. 

In Fig. 8 we show the magnetizations for the NFL states of various parameter cases. It is 
normalized by a quantity, M*, where M* is the magnetization under very weak magnetic field at 
a temperature, T* . The quantity, T* is defined as the temperature that — InT dependence of the 
magnetization begins. The strength of the — InT term is approximately characterized by 1/T*. 
The symbols shows the magnetizations for the f 2 - f 3 dominant fluctuation case, and the lines for 
the f 2 - f 1 dominant fluctuation one. As seen from the figure the normalized magnetizations under 
the same normalized magnetic field, H z /T* have almost the same temperature dependence. When 
the magnetic field is very weak, the temperature region that the magnetization has the — In T 
dependence spreads from T* to H z /10. The model of the f 2 - f 1 dominant fluctuation also gives 
sizable — InT divergence. 

§5. Magneto Resistance 

In the previous sections we have shown that the low energy spectrum of the present model is 
formed by the combination of the NFL and the LFL parts. We note that the phase shift which is 
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Fig. 6. Temperature dependence of the normalized magnetizations (a) and the differential of them by the logarithm 
of the temperature (b) for several magnetic fields, H z = 1 x 10~ 9 (o, solid line), 2 x 10 _9 (o, dashed line), 4 x 
10 -9 (A, dot-dashed line), 1 x 10 -8 (n, two-dots-dashed line), and 1 x 10~ 7 (+, three-dots-dashed line), respectively. 
The parameters are chosen so that the valence fluctuation between the f 2 and f 3 configurations becomes dominant: 
e(V^ ] ) = -1.2, e(r 6 ) = -1.05, e(4 2) ) = -0.8, U = 0.3, I = 8 and V 2 /2 = 0.035. The occupation number of the 
f-electron for H z = 1 x 10 -9 is 2.23. A = 8 is used, and about 400 states are retained at each renormalization step. 



calculated from the spectra of LFL part becomes very small at low energy fixed point. Therefore 
the electric resistivity caused by the LFL part will decrease with decreasing temperature. In this 
section we study the electric resistivity due to the NFL part. To simplify the analysis we use the 
fictitious model which is introduced in §3.2. When the magnetic field is applied, the spectrum at 
the low energy fixed point follows the LFL theory. The single particle excitations are given by the 
following effective Hamiltonian, 

H eS = H° L + Y,tf (r7)/rVr 7 
r 7 

+ ]T y/2Tf (r 7 ) /vr (fl lS0 r 7 + h.c.) , (5.1) 

r 7 
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Fig. 7. FCEL under the magnetic field, H z — 1 x 10 9 for the odd renormalization steps L. The parameters are 
the same as those in Fig. 6. See the caption of Fig. 1. 
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Fig. 8. Temperature dependence of scaled magnetizations for several magnetic fields. A = 8 is used, and about 400 
states are kept in each renormalization step. The parameters for the magnetizations, which are indicated by the 
symbols, are e(r^ 1) ) = -1.2, e(r 6 ) = -1.05, e(rf ) ) = -0.8, U = 0.3, 1 — 8 and V 2 /2 = 0.035, and H z /T* = 
1.35xl0" 2 (o), 2.7xl0 _2 (o), 5.4xlO~ 2 (A) and 0.135(D), where T* = 7.42 x 10" 8 . The occupation number of the f- 
electron for H z /T* = 1.35 x 10" 2 is 2.23. The magnetizations are normalized by M* = M(T = T* , H z = 1 x 10~ 9 ) 
= 1.07 x 10~ 3 . The parameters for the magnetizations which are indicated by the lines are ^F^ 1 ') = —0.9, 
e(r 6 ) = -0.75, e(rf ) ) = -0.5, U = 0.6, / = 8 and V 2 /2 = 0.035, and H z /T* = 5.5 x 10 -3 (solid line), 
1.1 x 10~ 2 (dashed line), 2.2 x 10~ 2 (dot-dashed line), 5.5 x 10~ 2 (two-dot-dashed line) and 0.55 (three-dot- 
dashed line), where T* = 1.82 x 10~ 8 . The occupation number of the f-electron for H z /T* = 5.5 x 10 -3 is 1.91. 
M* = M(T = T*,H Z = 1 x 10" 10 ) = 4.74 x 10" 3 . 
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where 5^ (r 7 ) and T e ^ (r 7 ) are the effective single particle level and effective hybridization width, 
respectively.© The effective quantities depend on the T and the suffix 7 specifying its component, 
and are determined by using the low energy excited states. The phase shift is given as, 



77^7) = - + arctan — ^4 — — . (5.2) 

/v ; 2 rf (r 7 )/A A v ; 




Figure 9 shows the magnetic fields dependence of the phase shift. The hybridization strength 
is increased from Figs, (a) to (c). We note even though the hybridization is large in Fig. (c), the 
low energy spectrum follows the NFL model for the case of no magnetic field. Each symbol shows 
the phase shift for each component of the T-irreducible representation. When the hybridization is 
weak, phase shift has the values about 7r/4 and 3tt/4 at weak magnetic field limit. This behavior 
is same to that of S'=l/2-TCKM@ When the magnetic field is increased, the phase shift changes 
gradually, and the strength of the scattering amplitude for all components decreases. In the case of 
the intermediate hybridization, the phase shift scarcely depends on the magnetic field. When the 
magnetic field is increased in the case of the large hybridization, the phase shift of three components 
changes to increase the strength of the scattering amplitude, while that of one component changes 
to decrease it. 

We calculate the magneto resistance by assuming that four components contribute with equal 
weight, 



Figure 10 shows the magneto resistance at T = 0, where the circle, diamond and triangle are 
calculated from the phase shift shown in Figs. 9(a), (b) and (c), respectively. Each resistance goes 
to 4.0 in the weak magnetic field limit. This value is half of the unitarity limit value. When the 
hybridization is small, we have the negative magneto resistance, while we have the positive magneto 
resistance for the large hybridization case. The magnetization for the case of large hybridization is 
shown in Fig. 11. We can see the — InT dependence with the large coefficient. The coefficient of the 
— In T term is large also for the small hybridization case. When the hybridization is intermediate, 
the resistance hardly depends on the magnetic field, and the temperature, T*, becomes high. This 
indicates that the initial; namely, the non-renormalized model is close to the low energy fixed point 
model. We note that the coefficient of the — In T term is small in this case. 

We note that the positive magneto resistance is recently observed at low temperature for 
U a Thi„ x Ru2Si20 For the positive magneto resistance case we can expect that the electric re- 
sistivity decreases with decreasing temperature. In fact, Affleck et al. have demonstrated that the 
magneto resistance is positive for the 5=l/2-TCKM with very strong exchange coupling case. 
In this case we can expect that the resistivity decreases with decreasing temperature. The renor- 
malized coupling decreases with decreasing temperature for the strong coupling model because the 
low energy fixed point of the NFL state is characterized by the intermediate coupling value. We 



i?(^) = ^(l-cos2r ? (r 7 )). 



(5.3) 



19 




.75 
.50 
.25 
.00 



VS ' 


(!) 1 


U) 1 




bl 


■ A 


A 


A 


A 




o 


❖ 


O 


O 




□ 


□ 


□ 


□ 

o 






10 -11 1Q 10 jQ-9 jQ-o JQ-, 1Q - 



Fig. 9. Magnetic field dependence of the phase shift. The parameters are e(r^) = —1.2, e(re) 

-0.8, U = 0.3, I = 8 and V 2 /2 = 0.035 (a), 0.13 (b) and 0.16 (c). The occupation number of the f-electron for 
H z = are 2.11(a), 2.27(b) and 2.29(c), respectively. Each symbol shows the phase shift for each component of 
the T-irreducible representation; [!>,+] (o), [r 7 , -] (o), [r 6 ,+] (A) and [r 6 , -] (□). 



-1.05, £ (rf J ) 



can also expect the positive magneto resistance when the low energy fixed point of the NFL state 
is approached from the stronger coupling side with renormalization step. The larger hybridization 
case of the present model seems to show this situation. We note that the present model can give 
the positive magneto resistance together with the sizable — In T term of the magnetic susceptibility. 
However, the S'=l/2-TCKM has very small coefficient of the — InT term with magnitude, 1/D, 
because the exchange constant must have comparable value with the band width to get the positive 
magneto resistance. 

From the calculations shown in this section we can expect that the anomalous properties of 
Ua;Thi_ x Ru2Si2 can be explained within the scenario of the NFL anomaly of the TCKM type. 
Hear we note the temperature dependence of the resistivity of the ETCAM. The resistivity of the 
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Fig. 10. Magneto resistance for the several cases of hybridization strength. The parameters are the same as those 
in Fig. 9 a (o), b (o) and c (A). 
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Fig. 11. Temperature dependence of the normalized magnetizations for the positive magneto resistance case. The 
parameters are the same as those in Fig. 9 (c). The magnetic fields are H z = 1 x 10~ 12 (o), 2 x 10~ n (o), 4 x 10~ 10 (A) 
and 1 x 10 _9 (n), respectively. 



ETCAM decreases with decreasing temperature, when the large hybridization is assumed.E-3' In 
the very low energy region of the NFL state, the scattering amplitude decreases to half of the 
unitarity limit value from the lager value (near the unitarity limit) of the preceding energy region. 
The preceding region is characterized by the Kondo effect of the two-channel Anderson model with 
larger hybridization. The behavior of the FCEL of the ETCAM have similar characteristics to 
that of the present model with larger hybridization. It shows the temperature dependence of the 
magnetization and the resistivity similar to that of U x Thi_ x Ru2Si2. The ETCAM seems to be 
an effective model of the NFL part of the realistic model at low temperature. The local spin in 
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the ETCAM will correspond to the f-electron with Iy symmetry, while the two channels will 

(2) 

correspond to electrons with Tj and Tq symmetries of the present model. However, at present it 
is not easy to relate the interaction terms between two models. 

§6. Stability of non-Fermi-liquid State 

In the above sections we have shown the NFL behavior of the system with the non-Kramers 
doublet CEF ground state. Because the low energy spectrum of the NFL part has the same form 
as that of S=l/2-TCKM, an entropy at the zero temperature is expected to remain. However, 
it conflicts with the law of thermodynamics. The residual entropy should be released by some 
remaining mechanisms. We have shown when the magnetic field is applied, the NFL state breaks 
below the low energy, H z /10, so the residual entropy is released. 
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Fig. 12. Temperature dependence of the magnetizations (a) and the differential of them by the logarithm of tem- 
perature (b) under the orthorhombic CEF. The magnetic field H z = 1 X 10 -10 is applied. The split of the two 
singlets of f 2 configuration is; 8 — (o, solid line), 2.2 x 10 -9 (o, dashed line), 2.1 x 10 -8 (A, dot-dashed line), 
2.1 x 10 -7 (□, two-dots-dashed line), respectively. The occupation number of the f-electron for 8 = is 1.91. 

We study an effect of the lowering of the crystal symmetry as another possible mechanism. When 
the crystal symmetry changes from the tetragonal type to the orthorhombic one, the non-Kramers 
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doublet CEF ground state splits into two singlets. If the split is large enough, the low energy 
states is expected to follow the LFL model, and the entropy at low temperature will go to zero. 
In Fig. 12 (a) we show the normalized magnetizations for several cases of the orthorhombic CEF. 
When the split, S, is zero, the result corresponds to that for the tetragonal CEF. Because the applied 
magnetic field, H z , is very weak, the magnetization shows — InT dependence for wide temperature 
region; from T* ~ 10~ 8 to H z /10 ~ 10~ n . The — InT dependence of the magnetic susceptibility is 
suppressed in the low temperature region, T < (5/10. But for the small split, the — InT dependence 
is still shown for the temperature region from T* to 5/10 as shown in Fig. 12 (b). The region is 
decreased with increasing 5, and it disappears for 5 ~ T*. The saturation of the magnetization 
occurs in the temperature range where the energy spectrum is explained by the LFL Hamiltonian. 
The residual entropy will be released at the temperature that the magnetization saturates. If the 
lowering of the crystal symmetry is not large, the NFL behavior will be observed in a restricted 
temperature region. 

§7. Summary 

We have investigated the NFL state of the IAM with non-Kramers doublet ground state of the f 2 
configuration under the tetragonal CEF. The low energy spectra, the temperature dependence of 
the magnetizations and the magneto resistances of the model are calculated by employing the NRG 
method. The spectra are explained by the combination of the NFL and the LFL parts which are 
independent with each other. One orbital of electrons with the symmetry contributes to the 
spectrum of LFL part, and the other two orbitals contribute to that of the NFL part. The phase 
shift of the conduction electron for the LFL component becomes very small at very low energy. 
The NFL part of the spectra has the same form to that of S=l/2-TCKM which has been derived 
by the CFT. 

The magnetization under the weak magnetic field shows the — InT dependence, and it saturates 
below the temperature, H z /10. The coefficient of the — InT term has magnitude of 1/T*, where 
T* is the temperature that the — InT dependence begins. The saturation is consistent with the 
break of the NFL state under the magnetic field. The magnetization is well scaled by T* and M* 
= M(T = T* , H z /T* <C 1). The magnetic field dependence of the phase shift is calculated from 
the low energy spectrum. When the hybridization is weak, the scattering amplitude decreases with 
increasing the magnetic field. But for large hybridization, the scattering amplitude increases. The 
existence of the positive magneto resistance suggests that the resistivity decreases with decreasing 
temperature. We have the parameter region for the IAM that the — InT term of the magnetic 
susceptibility has the sizable coefficient, and the resistivity decreases with decreasing temperature. 
There is the possibility that the anomalous properties of U x Thi_ x Ru2Si2 can be explained by the 
NFL scenario of the TCKM type for the present IAM. However, Actual U ion obeys L-S coupling 
scheme, so we need to consider complicated structures of the electronic states of the multi-electron 
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configurations for the quantitative discussions. 

The lowering of the crystal symmetry breaks the NFL behavior at the temperature about (5/10, 
where 5 is the orthorhombic CEF splitting. The residual entropy is released at around this tem- 
perature. If 5 is small, the NFL behavior is still expected in the intermediate temperature range. 

Acknowledgments 

We would like to thank H. Amitsuka, M. Koga, H. Shiba and Y. Kuramoto for useful discussions 
and comments. This work has been partly supported by Grant-in-Aid Nos. 06244104, 09640451 
and 09244202 from the Ministry of Education, Science and Culture of Japan. The computation has 
been done at 'the Computer Center, Tohoku University', 'the facilities of the Supercomputer Center, 
Institute for Solid State Physics, University of Tokyo' and 'the Computer Center of Institutes for 
Molecular Science, Okazaki National Research Institutes'. 



H. Amitsuka, T. Hidano, T. Honma, H. Mitamura and T. Sakakibara: Physica B 186-188 (1993) 337. 
H. Amitsuka and T. Sakakibara: J. Phys. Soc. Jpn. 63 (1994) 736. 

H. Amitsuka, T. Sakakibara, A. de Visser, F. E. Kayzel and J. J. M. France: Physica B 230-232 (1997) 613. 
D. L. Cox: preprint, |eond-mat /9704103 , to be published in Adv. Phys. 
D. L. Cox: Phys. Rev. Lett. 59 (1987) 1240. 
D. L. Cox: Physica B 186-188 (1993) 312. 
P. Nozieres and A. Blandin: J. Phys. (Paris) 41 (1980) 193. 
A. M. Tsvelick and P. B. Wiegmann: J. Stat. Phys. 38 (1985) 125. 
P. D. Sacramento and P. Schlottmann: Phys. Rev. B 43 (1990) 295. 

I. Affleck and A. W. W. Ludwig: Phys. Rev. B 48 (1993) 7297. 
O. Sakai, Y. Shimizu and N. Kaneko: Physica B 186-188 (1993) 322. 
H. Kusunose, K. Miyake, Y. Shimizu and O. Sakai: Phys. Rev. Lett. 76 (1996) 271. 
O. Sakai, S. Suzuki, Y. Shimizu, H. Kusunose and K. Miyake: Solid State Commun. 99 (1996) 461. 
S. Suzuki, O. Sakai and Y. Shimizu: Solid State Commun. 104 (1997) 429. 
M. Koga and H. Shiba: J. Phys. Soc. Jpn. 64 (1995) 4345. 
M. Koga and H. Shiba: J. Phys. Soc. Jpn. 65 (1996) 3007. 
M. Koga and H. Shiba: J. Phys. Soc. Jpn. 66 (1997) 1485. 

0. Sakai, S. Suzuki and Y. Shimizu: Solid State Commun. 101 (1997) 791. 

Y. Shimizu, O. Sakai and S. Suzuki: The analysis of the low energy spectrum will be partly published in J. Mag. 
& Mag. Mat. 

P. H. Butler: Point Group Symmetry Applications, Methods and Tables (Plenum Press, New York, 1981). 
H. R. Krishna-Murthy, J. W. Wilkins and K. G. Wilson: Phys. Rev. B 21 (1980) 1003. 

1. Affleck, A. W. W. Ludwig, H. B. Pang and D. L. Cox: Phys. Rev. B 45 (1992) 7918. 

The small difference may be due to the residual small interaction between the NFL and the LFL parts because 
we listed the energies at L — 31 step. 
[24] H. Amitsuka: private communication. 



24 



